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A phenomenological momentum-independent (MID) model is constructed to describe the equation 
of state (EOS) for isospin asymmetric nuclear matter, especially the density dependence of the 
nuclear symmetry energy -Bsym (/!>)• This model can reasonably describe the general properties of the 
EOS for symmetric nuclear matter and the symmetry energy predicted by both the sophisticated 
isospin and momentum dependent MDI model and the Skyrme-Hartree-Fock approach. We find 
that there exists a nicely linear correlation between i^sym and L as well as between Jo/Ko and Ko, 
where L and Ksym represent, respectively, the slope and curvature parameters of the symmetry 
energy at the normal nuclear density po while Ko and Jo are, respectively, the incompressibility 
and the third-order derivative parameter of symmetric nuclear matter at po- These correlations 
together with the empirical constraints on Kq, L and iJsym(po) lead to an estimation of —477 
MeV < ifsat,2 < —241 MeV for the second-order isospin asymmetry expansion coefficient for the 
incompressibility of asymmetric nuclear matter at the saturation point. 

PACS numbers: Equation of state of nuclear matter, isospin, the symmetry energy 



I. INTRODUCTION 

The study of the isospin degree of freedom in nuclear 
physics has recently attracted much attention due to the 
establishment of many radioactive beam facilities around 
the world. Besides the many existing radioactive beam 
facilities and their upgrades, such as the Cooling Storage 
Ring (CSR) facility at HIRFL in China d, many more 
are being constructed or under planning, including the 
Radioactive Ion Beam (RIB) Factory at RIKEN in Japan 
0, the FAIR/GSI in Germany [|, SPIRAL2/GANIL 
in France [3], and the Facility for Rare Isotope Beams 
(FRIB) in the USA These new facilities offer the 
possibility to study the properties of nuclear matter or 
nuclei under the extreme condition of large isospin asym- 
metry. The ultimate goal of such study is to extract 
information on the isospin dependence of in-medium nu- 
clear effective interactions as well as the equation of state 
(EOS) of isospin asymmetric nuclear matter, particularly 
its isospin-dependent term or the density dependence of 
the nuclear symmetry energy. This knowledge, especially 
the latter, is important for understanding not only the 
structure of radioactive nuclei, the reaction dynamics in- 
duced by rare isotopes, and the liquid-gas phase transi- 
tion in asymmetric nuclear matter, but also many critical 
issues in astrophysics [1, 0, i, i, [lO, EH, E [H, Q [11 ■ 

The EOS of nuclear matter is one of fundamental ques- 
tions in nuclear physics. For symmetric nuclear matter, 
the EOS is relatively well-determined after about more 
than 30 years of studies in the nuclear physics community. 
The incompressibility of symmetric nuclear matter at its 
saturation density po has been determined to be 240 ± 20 
MeV from the nuclear giant monopole resonances (GMR) 
EH EE, [13 and the EOS at densities of 2po < p < 
5po has also been constrained by measurements of col- 
lective flows in nucleus-nucleus collisions [1] and of sub- 
threshold kaon production [2l|,[2^ in relativistic nucleus- 



nucleus collisions. On the other hand, for asymmetric nu- 
clear matter, the EOS, especially the density dependence 
of the nuclear symmetry energy, is largely unknown. Al- 
though the nuclear symmetry energy at po is known to 
be around 30 MeV from the empirical liquid-drop mass 
formula 23 , [13] , its values at other densities are poorly 
known [fi 1? j . Various microscopic and phenomenological 
models, such as the relativistic Dirac-Brueckner-Hartree- 
Fock (DBHF) [H, Hi, [13, [H, [H, [10, [m and the non- 
relativistic Brueckner-Hartree-Fock (BHF) [H, [H, Hi] 
approach, the relativistic mean-field (RMF) model based 
on nucleon-meson interactions (l^ . [35 , 36, 37], and the 
non-relativistic mean-field model based on Skyrme-like 
interactions [H, [M [10, [ill, [13, [H, [13, [H , have been used 
to study the isospin-dependent properties of asymmetric 
nuclear matter, such as the nuclear symmetry energy, the 
nuclear symmetry potential, the isospin-splitting of nu- 
cleon effective mass, etc., but the predicted results vary 
widely. In fact, even the s ign of the symmetry energy 
above 3po is uncertain [46l . w\ . The theoretical uncer- 
tainties are mainly due to the lack of knowledge about the 
isospin dependence of in-medium nuclear effective inter- 
actions and the limitations in the techniques for solving 
the nuclear many-body problem. 

In the present work, we construct a phenomenologi- 
cal momentum-independent (MID) model which can rea- 
sonably describe the general properties of symmetric nu- 
clear matter and the symmetry energy predicted by both 
the sophisticated isospin and momentum dependent MDI 
model and the Skyrme-Hartree-Fock approach with dif- 
ferent Skyrme forces. In particular, the density func- 
tional of the symmetry energy constructed in the MID 
model is shown to be very flexible and can mimic very 
different density behaviors by varying only one param- 
eter. We find that there exists a nicely linear correla- 
tion between i^sym and L as well as between Jq/Kq and 
Kq, where L and Ksym represent, respectively, the slope 
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and curvature parameters of the symmetry energy at the 
normal nuclear density po while Kq and Jq are, respec- 
tively, the incompressibility and the third-order deriva- 
tive parameter of symmetric nuclear matter at po- These 
correlations together with the empirical constraints on 
Kq, L and £'sym(po) lead to an estimate of —477 MeV 
< Ksiit.2 < —241 MeV for the second-order isospin 
asymmetry expansion coefficient for the incompressibil- 
ity of asymmetric nuclear matter at the saturation point, 
which is presentl y l argely uncertain and being heavily 
discussed [H, M, M El, HO] ■ 

The paper is organized as follows. In Section [Til we 
discuss general properties of asymmetric nuclear mat- 
ter, and then introduce the momentum independent MID 
model. The results and discussions are presented in Sec- 
tion mil A summary is then given in Section IIVI 



II. THEORETICAL MODELS AND METHODS 

A. Equation of state of asymmetric nuclear matter 

The EOS of isospin asymmetric nuclear matter, given 
by its binding energy per nucleon, can be expanded to 
2nd-order in isospin asymmetry S as 



Eip, S) = Eoip) + Esy^ip)d^ + 0{5^ 



(1) 



where p = p„ -I- pp is the baryon density with p„ and 
Pp denoting the neutron and proton densities, respec- 
tively; 5 = {pn — Pp)/{Pp + Pn) is the isospin asymmetry; 
Eq{p) ~ E{p,S ~ 0) is the binding energy per nucleon 
in symmetric nuclear matter, and the nuclear symmetry 
energy is expressed as 



1 d^E{p,S) 
2! 86^ 



\s=o- 



(2) 



The absence of odd-order terms in d in Eq. ([T]) is due to 
the exchange symmetry between protons and neutrons in 
nuclear matter when one neglects the Coulomb interac- 
tion and assumes the charge symmetry of nuclear forces. 
The higher-order coefficients in S are usually very small 
and negligible, e.g., the magnitude of the S'^ term at nor- 
mal nuclear density po is estimated to be less than 1 MeV 
in microscopic many-body approaches [5l|, [13, HI] ■ Ne- 
glecting the contribution from higher-order terms in Eq. 
(fT|) leads to the well-known empirical parabolic law for 
the EOS of asymmetric nuclear matter, which has been 
verified by all many-body theories to date, at least for 
densities up to moderate values . As a good approxi- 
mation, the density-dependent symmetry energy Egymip) 
can thus be extracted from the parabolic approximation 
of E,y„,{p) « E{p, S=l)- E{p, S = 0). 

Around the normal nuclear density po, the binding en- 
ergy per nucleon in symmetric nuclear matter Eq [p) can 
be expanded, e.g., up to 3rd-order in density as 



Eo{p)=Eq{Po) 



Ko 
2! 



r + §x' + 0(x^), (3) 



where x is 9- dimensionless variable characterizing the 
deviations of the density from the saturation density pQ 
of the symmetric nuclear matter and it is convention- 
ally defined a,s x = [p ~ Po)/3/Oo- Eq{pq) is the binding 
energy per nucleon in symmetric nuclear matter at the 
saturation density po the other coefficients can be 
calculated as 



2<PEo{p) 



dp^ 



\p=pg,J0 - Z/Po 



dp^ 



(4) 



Obviously, there is no linear x term in Eq. ([3]) according 
to the definition of the saturation density po- Kq is the 
incompressibility coefficient of symmetric nuclear matter 
and it characterizes the curvature of Eq{p) at po. The 
coefficient Jo corresponds to the third-order derivative 
parameter of symmetric nuclear matter at po. In the lit- 
erature, people usually neglect the higher-order terms in 
Eq. ([3]) and obtain the following parabolic approximation 
to the EOS of symmetric nuclear matter: 



Eoip) = ^o(po) 



'^x' + oix'). 



(5) 



Similarly, around po, the nuclear symmetry energy 
Esyniip) can be expanded, e.g., up to 2nd-order in density 
as 



Esymip) = Esymipo) + Lx + 



sym 2 
A 



2! 



where L and Ksym are the slope parameter and curvature 
parameter of the nuclear symmetry energy at po, i.e.. 



3po — i— 

dp 



(P) 



P=Po ' -^sym 



2 d'^Esymjp) I 
dp^ ' 



(6) 

The coefficients L and A'sym characterize the density de- 
pendence of the nuclear symmetry energy around normal 
nuclear density po, and thus carry important information 
on the properties of nuclear symmetry energy at both 
high and low densities. 

The incompressibility is an essential quantity for nu- 
clear matter and conventionally it is defined at the satu- 
ration density psat for asymmetric nuclear matter where 
we have P{p, S) ^ (the incompressibility coefficient at 
the saturation density is called isobaric incompressibility 
coefficient in [s^]) and thus it can be expressed as 



K,,ti6) = 9pl 



P=Ps: 



(7) 



For asymmetric nuclear matter, the isobaric incompress- 
ibility coefficient Kgnti^) can be expressed up to 2nd- 
order in d as ISa 



with 



Arsat,2 = A'sym — 6L — -r^L. 



Ko 



(8) 



(9) 
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The coefficient i^sat.2 essentially reflects the isospin de- 
pendence of the isobaric incompressibility of asymmetric 
nuclear matter. 

If we use the parabolic approximation to EOS of sym- 
metric nuclear matter, i.e., Eq. (O, then the i^sat,2 is 
reduced to 



6L 



(10) 



and this expression has been extensively used to charac- 
terize the isospin dependence of the incompressibility of 
asy mmetric nuclear matter in the literature [56l . [stI [ssI . 
|59|. Obviously, we have 



-K^sat,2 = Kasy —L, 



(11) 



and thus the coefhcient K^sy could be a good approxima- 
tion to A'sat,2 if ^0 is negligible or the slope parameter of 
the symmetry energy L is very small. 

It is believed that information on -R'sat,2 can in princi- 
ple be extracted experimentally by measuring the GMR 
in neutron-rich nuclei [ssj . Usually, one can define a finite 
nucleus incompressibility Ka{N, Z) for a nucleus with N 
neutrons and Z protons {A = N + Z) hy the energy of 
GMR^Jgmr, i.e., 



GMR 



Ih^KAjN, Z) 



(12) 



where m is the nucleon mass and ^r^^ is the mean square 
mass radius of the nucleus at ground state. Similar to 
the semi-empirical mass formula, the finite nucleus in- 
compressibility Ka{N, Z) can be expanded as 



N - Z 
A 



-Kco\i\- 



(13) 

where Kq, ii'surf, Kt, and i^coui represent the volume, 
surface, symmetry, and Coulomb terms, respectively. 
The Kr parameter is usually thought to be equivalent 
to the Ksat,2 parameter. It should be noted here that 
the -ftrsat,2 parameter is theoretically a well-defined phys- 
ical quantity while the value of the K^- parameter may 
depend on the detailed truncations in the expansion sim- 
ilarly to the semi-empirical mass formula. Earlier at- 
tempts based on the above method have given widely 
different values for the Kr parameter. For example, a 
value of K-T = —320 ± 180 MeV with a large uncertainty 
was obtained in Ref. [60j from a systematic study of 
the GMR in the isotopic chains of Sn and Sm. In this 
analysis, the value of Kq was found to be 300 ± 25 MeV, 
which is somewhat larger than the commonly accepted 
value of 240 ± 20 MeV. In a later study, an even less 
stringent constraint of —566 ± 1350 < Kt < 139 ± 1617 
MeV was extracted from the GMR of finite nuclei, de- 
pending on the mass region of nuclei and the number of 
parameters used in parameterizing the incompressibility 



of finite nuclei [6l| . Most recently, a much stringent con- 
straint of Kt = —550 ± 100 MeV has been obtained in 
Ref. [H, [III from measurements of the isotopic depen- 
dence of the GMR in even-A Sn isotopes. 



A phenomenological momentum-independent 
MID model 



In the present work, we will mainly use three 
models, i.e., the isospin and momentum depen- 
dent MDI interaction [3^, the Hartree-Fock approach 
based on Skyrme interactions, and a phenomenological 
momentum- independent interaction (MID). The MDI in- 
teraction [3^ is based on the finite-range Gogny effective 
interaction and has been used extensively in the litera- 
ture [iBl ■ The SHF approach is a well-known mean-field 
theory and has been extensively used in the literature 
for its simplicity. A very useful feature of these models is 
that analytical expressions for many interesting physical 
quantities in asymmetric nuclear matter at zero temper- 
ature can be obtained. Here we only introduce the MID 
model and for the MDI and SHF models, one can refer 
to, e.g., Refs. [38l.[6^. 

In the momentum-independent MID model, following 
the results from SHF approach with the zero-range and 
momentum-independent Skyrme interaction, the poten- 
tial energy density Vmid(p, (5) of a cold symmetric nu- 
clear matter at total density p and isospin asymmetry 6 
is parametrized as 



a p'^ (3 p'^+i 
Vmid(p,5) = -z— + — — — — 

2 po cr + 1 po'^ 



+ pEP^"Up)S\ (14) 



In the MID model, the 4th-order and higher-order nu- 
clear symmetry energy are not included and we assume 
they can be negligible. The parameters a, (3 and a are 
determined by the binding energy per nucleon Eq{pq) = 
— 16 MeV and the incompressibility Kq at the saturation 
density po = 0.16 fm"'^ 



a = -29.47-46.74 
P ^ 23.37 



Kn + 44.21 



Ko - 166.11 
Ko + 254.53 



Ko - 166.11 
A'o + 44.21 
210.32 ' 



(MeV), (15) 
(MeV), (16) 
(17) 



where the unit of Kq is MeV 
in the MID model, it is parametrized as 



For the potential part of the symmetry energy i?|'°4(p) 



EP^^^ip) = i?r;*n(po)(i - y)- + yE!;M , 



with EP°l^ipo 



Esyniipa) - E^™{po) 



(18) 
17.7 MeV 



following Esymipo) = 30 MeV and E^^^ipo) = 
& (^Po)^^^ = 12.3 MeV. The default value of the 
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7sym parameter is taken to be 4/3 in the MID model 
following the -Esym(/f) in the MDI interaction (we will 
see how the 7syni parameter affects the symmetry energy 
in the following). Similarly to the x parameter intro- 
duced in the MDI interaction jsl], the dimensionless y 
parameter is introduced to mimic various Esym{p) pre- 
dicted by different microscopic and/or phcnomcnological 
many-body theories for a fixed 7sym parameter. As we 
will show later, for 7syni 4/3, adjusting the y value 
can nicely reproduce the -Esym(p) in the MDI interaction 
with X = —1, 0, and 1. 

In the MID model, the EOS of symmetric nuclear mat- 
ter can thus be written as 



100 



Eoif) = 



3^2 /'Stt 



2/3 



a p 
2^ 



/3 



lOm V 2 / 2/90 0- + 1 
and the symmetry energy can be expressed as 



(19) 



2/3 



+ [£;.y„,(po)-i?syS.(po)](l-2/)^ 



Po 



Po 



^symvA^u; -^symV 

which leads to 



(20) 



pkin . 
^sym \ 

-3y(7 



sym - 1) [Esyniipa) " E^^^^ipo)] (21) 



sym 



syii 

ki' 

sym V 
ki 

symV 



92/7sym(7sym - 1) [Esymipo] ~ -E^ym(Po)] 
-2i?s'ym(P0). (22) 



III. RESULTS AND DISCUSSIONS 



The nuclear symmetry energy and correlation 
between L and Jfsvm 



As mentioned above, in the MID model, the density 
dependence of the symmetry energy can be adjusted by 
varying the y parameter as shown in Eq. (|20p . As an 
example, we show in Figure [1] the density dependence 
of the symmetry energy from the MID interaction with 
y = —3.4, —0.73, and 1.8. The corresponding results 
from the MDI interaction with x = 1, 0, and —1 as well as 
the widely used APR (Akmal-Pandharipandc-Ravcnhall) 
prediction (63j arc also included for comparison. Indeed, 
one can see that the MID interaction can give a nice 
description on the density dependence of the symmetry 
energy predicted by the sophisticated MDI interaction 
from the very soft {x = 1) to the very stiff one {x — I). 
Furthermore, it is seen that the APR prediction for the 
symmetry energy at subsaturation densities lies right be- 
tween that with X — and —1, and especially the sym- 
metry energy with x — (and y ~ —0.73) resembles very 



~^ \ ' \ ' r 

MDI(x): Thick lines , 
IVIID(y): Jh\r\ lines .' 




FIG. 1: (Color online) Density dependence of the symmetry 
energy from the phenomenological MID interaction predic- 
tion with y — —3.4, —0.73, and 1.8. The results from the 
MDI interaction with a:: = 1, 0, and —1 and the widely used 
APR (Akmal-Pandharipande-Ravenhall) prediction [63| are 
also included for comparison. 



well the APR prediction up to about 3.5po- These fea- 
tures imply that the density functional of the symmetry 
energy shown in Eq. (j20p is very flexible and can give 
a quite general description for density dependence of the 
symmetry energy. 

The parameters L and -K^sym are determined by the 
density dependence of the symmetry energy around po. 
In recent years, significant progress has been made both 
experimentally and theoretically in extracting the infor- 
mation on the symmetry energy at sub-saturation den- 
sity from heavy-ion reactions. Using the isospin and 
momentum-dependent IBUU04 transport model with in- 
medium NN cross sections, the isospin diffusion data 
were found to be consistent with the symmetry energy 
from the MDI interaction with x between and — 1, 
which can be parametrized by Esymip) ~ 31.6(p/po)'*^ 
with 7 = 0.69 — 1.05 at subnormal density (p < po) 
[45l . [ssl . [6^ . [65| . and has led to the extraction of 61 
MeV < L < 111 MeV and -82 MeV < Ksym < 101 



MeV |45|,l58|,l6J,l65l. Using the Skyrme interactions con- 
sistent with the EOS obtained from the MDI interac- 
tion with X between and —1, the neutron-skin thick- 
ness of heavy nuclei calculated within the Hartree-Fock 
approach is consistent with available experimental data 
[45I . [66j and also that from a rclativistic mean-field model 
based on an accurately calibrated parameter set that 
reproduces the GMR in ^"Zr and ^°*Pb as well as the 
isovector giant dipole resonance of ^°^Pb [63|- The ex- 
tracted symmetry energy further agrees with the sym- 
metry energy i?sym(p) = 31.6(p/po)°'^^ recently obtained 
from the isoscaling analyses of isotope ratios in interme- 
diate energy heavy ion collisions [68j , which gives L w 65 



MeV and Ks- 



-61 MeV. Furthermore, it is interest- 



ing to mention that the above limited range of i?sym(p) 
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at subsaturation density is essentially consistent with the 
symmetry energy E^y^{p) = 12.5(p/po)^/^ + n.Q{plpoV 
with 7 = 0.4 — 1.05, extracted very recently from analy- 
ses using the ImQMD (Improved QMD) model which can 
reproduce both the isospin dilfusion data and the double 
neutron/proton ratio simultaneously [6^ . The symme- 
try energy E,y^{p) = 12.5(p/po)'/^ + 17.6(p/po)'' with 
7 = 0.4 — 1.05 thus leads to the constraints of 46 MeV 
< L < 80 McV and -82 MeV < A'^ym < -36 McV. 
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FIG. 2: (Color online) Correlation between i^sym and L from 
the MID interaction with 7sym = 4/3 and 5/3, the MDI in- 
teraction and the SHF prediction with 63 popular Skyrme 
forces. 

It should be noted that all the above constraints on 
L and Kgy-^ are based on some unique energy density 
functionals and thus special correlation between L and 
Ksym has been implicitly assumed. It is thus interesting 
to see if there exists a universal correlation between L 
and Ksym- For the MDI interaction, the L and Kgym 
both change linearly with the parameter x and therefore 
they are linearly correlated by varying the parameter x 
[S^ . [70j . Similarly, for the MID interaction, one can see 
from Eq. dH]) and Eq. that the L and A'sym both 

change linearly with the parameter y, and thus they are 
also linearly correlated by varying the parameter y. In 
particular, we have 



ATsym = 37symi -|- E'^y™ (po) (37sym - 2) - 97sy„i£;sy„i (po) • 

(23) 

Also the L and ATsym are expected to be correlated within 
the SHF energy density functional. Shown in Figure [2] 
are the correlation between ATsym and L from the MID 
interaction with 7syin = 4/3 and 5/3 (E'gy™ (po) = 12.3 
MeV and £'sym(po) = 30 MeV), the MDI interaction and 
the SHE prediction with 63 popular Skyrme forces. The 
63 Skyrme forces include the 51 forces used in Ref. [t^I 
and 12 new forces, i.e., Z, E,,, E, Z;, SkSC4, SI, 
SII, SIII, SIV, SV, and SVI. All these Skyrme forces 
predict the saturation density and the symmetry en- 
ergy satisfying 0.140 fm""^ < po < 0.165 fm""^ and 25 
McV< Esym(po) < 37 McV, respectively. 



It is interesting to see that the A'sym parameter indeed 
displays approximately a linear correlation with the L 
parameter for the SHF prediction with the 63 Skyrme 
forces and this linear correlation is nicely reproduced 
by the MDI interaction and the MID interaction with 
7sym = 4/3. For the MID interaction, one can see from 
Eq. (1^^ that the 7sym parameter controls the shape 
(slope) of the linear correlation between L and -fCsym- 
Furthermore, it is seen from Figure [2] that there are a 
few Skyrme forces deviate from the linear correlation ob- 
tained by the MDI interaction and the MID interaction 
with 7syin = 4/3. In order to consider the uncertainty 
of the shape (slope) for the correlation between L and 
ATsym, we thus include the result with 7sym = 5/3 for 
the MID interaction. The correlation between Ksym and 
L from the SHF prediction with the 63 Skyrme forces 
is nicely consistent with that from the MID interaction 
with 7syin = 4/3 and 5/3. The linear correlation between 
-f^sym and L implies that one can obtain ATgym from L. 



B. Correlation between Jo and Kq 
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FIG. 3: (Color online) Jo and Jq/Kq as a function of Ko 
from the MID interaction, the MDI interaction and the SHF 
prediction with 63 popular Skyrme forces. 

While Kq has been relatively well determined, the Jo 
parameter is poorly known and actually there is no any 
experimental information on the Jq parameter. In the 
MID model, from Eq. one can easily calculate the 
Jo parameter as 



Jo = 27p; 



1 



70.1 



(Kl - 332.2^0 



4243.2) (MeV), (24) 



where the unit of Kq is MeV. Therefore, in the MID 
model, the Jq parameter is quadratically correlated with 
Kq. Shown in Figure [3] are Jo and Jq/Kq as func- 
tions of Kq. Also included in Figure [3] are the corre- 
sponding results from the MDI interaction and the SHF 
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prediction with the 63 Skyrme forces. It is interest- 
ing to see that the correlation between Jq and Kq is 
quite consistent for the three different models, namely, 
the MID interaction, the MDI interaction and the 63 
Skyrme forces in the SHF approach. In particular, the 
Jo/Kq displays approximately a linear correlation with 
Kq. This linear correlation can be easily understood from 
Eq. ([23]). On the r.h.s of Eq. (gll), the last term is 
very small compared with the first term and the second 
term and thus one has Jq « (^Kq — 332.2ifo)i and 
then Jq/Kq « {Kq — 332.2) with the unit of Kq be- 
ing MeV. We note here that the correlation between Jo 
and Kq obtained in the present work is also consistent 
with the early finding by Pearson [7l| . While there is no 
any empirical constraint on the Jo parameter, we assume 
in the present work the correlation between Jq and Kq 
from the MID interaction is valid and then we can obtain 
Jq/Kq from the experimental constraint on Kq. 

C. Phenomenological MID model constraint on the 
A'sat,2 parameter 
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FIG. 4: (Color online) A'aat,2 as a function of L from the 
MID interaction with 7sym = 4/3 (a) and 5/3 (b) for different 
values of Kq and -Bsym(po). The shaded region indicates con- 
straints within MID interaction with 220 MeV < Ka < 260 
MeV, 25 MeV < £'sym(po) < 35 MeV, and 46 MeV < L < 111 
MeV limited by the heavy-ion collision data. 

As shown in Eq. ([§]), the /irsat,2 parameter is com- 
pletely determined by Jq/Kq, L, and -K'sym- Based on 
the correlations shown in Figure [5] and Figure [31 we can 
now extract information on the Ksat,2 parameter from 
the experimental constraints on the Kq parameter and 
the L parameter within the MID model. As pointed out 
previously, the value of Kq has been relatively well de- 
termined to be 240 ± 20 MeV from the nuclear GMR 
[H, [13, [H, [H, [131. The slope parameter L has been 
found to correlate linearly with the neutron-skin thick- 
ness of heavy nuclei and thus can in principle be deter- 
mined from measured thickness of the neutron skin of 
such nuclei [H, [H, [zl, [zl, [zl, [ll, [ll, [t^ . Unfortunately, 
because of the large uncertainties in the experimental 
measurements, this has not yet been possible so far. The 



proposed experiment of parity- violating electron scatter- 
ing from ^"^Pb, i.e.. Parity Radius Experiment (PREx) 
at the Jefferson Laboratory is expected to give an inde- 
pendent and accurate measurement of its neutron skin 
thickness (within 0.05 fm) [tI, [?§]. On the other hand, 
as mentioned previously, heavy-ion collisions, especially 
those induced by neutron-rich nuclei, provide a unique 
tool to explore the density dependence of the symmetry 
energy and thus the L parameter. 

In the MID model, from Eqs. ^ and we have 

Ks.t,2 = -(^+6-37sym)L+(37sym-2)£:,^™(po) 

Ao 

-97sy,n£'sym(Po)- (25) 

Shown in Figure [4] is Ksiit,2 as a function of L from the 
MID interaction with 7sym = 4/3 (panel (a)) and 5/3 
(panel (b)) for Kq = 220, 240, and 260 MeV. In Eq. 
(PS)) . Jo can be obtained from Eq. pi)) for a fixed ATo 
value. In the MID interaction, E^^^{po) = 12.3 MeV and 
Asym(po) = 30 MeV have been used as a default. From 
Eq. (|23|1. one can see that the correlation of ATsym and L 
also depends on Esym{po)- To consider the uncertainty 
due to the £'sym(po)j we thus also include in Figure [4] the 
resuhs with A'o = 220 MeV and Esym{po) = 25 MeV as 
weh as ATo = 260 MeV and E^ymipo) = 35 MeV, which 
represent, respectively, the upper and lower boundaries 
for a fixed L. The shaded region in Figure [4] further con- 
siders the constrained L values from heavy-ion collision 
data, namely, 46 MeV < A < 111 MeV. The lower hmit of 
A = 46 MeV is obtained from the lower boundary of the 
ImQMD analyses on the isospin diffusion data and the 
double neutron/proton ratio [69| while the upper limit of 
A = 111 MeV corresponds to the upper boundary of A 
from the IBUU04 transport model analysis on the isospin 
diffusion data [H, [Bs, 64, 65i]. The constraint 46 MeV 
< A < 111 MeV is also consistent with the analyses of 
the pygmy dipole resonances [s^l, the giant dipole reso- 
nance (GDR) of ^°*Pb analyzed with Skyrme forces [Slj, 
the Thomas-Fermi model fitted very precisely to binding 
energies of 1654 nuclei [s^ . and the recent neutron-skin 
analysis [s^. These empirically extracted values for A 
represent the best and most stringent phenomenological 
constraints available so far on the nuclear symmetry en- 
ergy at sub-saturation densities. 

It is seen from Figure [3] that the Ksat,2 decreases with 
increasing A for 7sym = 4/3 while it increases with in- 
creasing A for 7sym = 5/3. This feature can be easily 
understood from Eq. ^t^. For 7syni = 4/3, Eq. (gSl) is 
reduced to 

Ks.t,2 = -(^ + 2)A - 12A,y„,(po) + 24.6 (MeV) (26) 
Ao 

while for 7sym = 5/3, it is reduced to 

A'sat,2 = -(^ + 1)A - 15A,y,„(po) + 36.9 (MeV). (27) 
Ao 
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For Kq = 240 ±20 McV, Jo/Kq can be found from Figure 
[3] (or Eq. dM])) to be from about -1.9 to -1.3. Therefore, 
Jsrsat,2 decreases (increases) with increasing L for 7sym = 
4/3 (5/3) following Eq. ^ (Eq. 

An interesting feature observed from Figure |4] is that 
the i^sat.2 parameter significantly depends on the sym- 
metry energy at the normal nuclear density i?sym(po)- 
This can be seen more clearly from Eqs. ([^5]) and ([77)1 
which indicate that changing -Esym(po) by 5 MeV leads 
to a variation of 60 — 75 MeV for Ksnt,2- This feature 
indicates that an accurate determination of -Esym(po) is 
important for determining the value of -firsat,2- From the 
shaded region indicated in Figure IH it is found that for 
7sym = 4/3, we have -429 MeV < Ksi,t.2 < -281 McV 
for L = 46 MeV while -477 MeV < ifsat,2 < -289 MeV 
for L = 111 MeV. For %yn, = 5/3, we have -476 MeV 

< A'sat.2 < -298 McV for L = 46 MeV while -459 McV 

< Arsat,2 < -241 McV for L = 111 MeV. These results in- 
dicate that within the MID model with the empirical con- 
straints of Ko = 240 ±20 MeV, 25 MeV < E,yrr,ipo) < 35 
MeV, and 46 MeV < L < 111 MeV, the A"sat,2 parameter 
can be varied from -477 MeV to -241 MeV. 

As shown in Eq. (|lip . the -fiTasy parameter corresponds 
to the ^Csat,2 parameter when Jq is zero, i.e., the parabolic 
approximation to the EOS of symmetric nuclear mat- 
ter Eq. ([5]) is valid. From the MID model, a vanish- 
ing Jo corresponds to a A'o value of about 340 MeV, 
which is significantly larger than the empirical value of 
240 ±20 MeV. In the MID model, we have Jo/Kn w -1.6 
for Ko = 240 MeV and thus Kss.t,2 ~ K^^y + 1-6L or 
Kasy ~ Arsat,2 " 1.6i. Therefore, the difference between 
Kiisy and Ksa,t,2 depends on L with a larger L value 
(stiffer symmetry energy) leading to larger difference. At 
this point, it should be stressed that the ifasy parame- 
ter is completely determined by the density dependence 
of the symmetry energy regardless of the EOS of sym- 
metric nuclear matter. Based on the IBUU04 transport 
model analysis on the isospin diffusion data [H, , a 
value of Ks,sy = —500 ± 50 MeV has been extracted 
from the symmetry energy obtained by the MDI inter- 
action with the X parameter between and —1. The 
constraint A'asy = —500 ± 50 MeV is quite consistent 
with the very recent constraint of K^sy ~ — 500^j^o MeV 
from the study of neutron skin of finite nuclei |83| . Fur- 
thermore, in the MDI interaction, we have —311 MeV 

< Kss^t,2 < -316 MeV from the prediction of the MDI in- 
teraction with the X parameter between and —1. There- 
fore, for the MDI interaction, the magnitude of ifsat,2 is 
significantly smaller than that of i^asy and is quite insen- 
sitive to the density dependence of the symmetry energy. 
These features indicate that the high-order Jq contribu- 
tion to Ksnt,2 generally cannot be neglected. 



IV. SUMMARY AND CONCLUSIONS 

We have constructed a phenomcnological momentum- 
independent MID model which can reasonably describe 



the general properties of symmetric nuclear matter and 
the symmetry energy predicted by both the sophisticated 
isospin and momentum dependent MDI model and the 
SHF approach with different Skyrme forces. In particu- 
lar, the density functional of the symmetry energy con- 
structed in the MID model is shown to be very flexible 
and can mimic very different density behaviors by vary- 
ing only one parameter. 



Based on the MID model, we have studied in detail 
the second-order isospin coefficient A'sat,2 which is deter- 
mined uniquely by L, Kgym and Jq/Kq. Our results in- 
dicate that the high-order Jq contribution to Arsat,2 gen- 
erally cannot be neglected, especially for larger L values. 
In addition, interestingly, it is found that there exists 
a nicely linear correlation between Kgym and L as well 
as between Jq/Kq and Kq for the three different mod- 
els used here, i.e., the MDI interaction, the MID inter- 
action, and the SHF approach with 63 Skyrme forces. 
From the MID model, the correlation between Ksym and 
L is further shown to depend significantly on the value 
of Esymipa)- These correlations and features enable us 
to extract the values of the Jo parameter and the ATsym 
parameter from the empirical information on Kq, L and 
Esymipo)- In particular, using the empirical constraints 
of A'o = 240 ± 20 MeV, 25 MeV < Asym(po) < 35 MeV, 
and 46 MeV < A < 111 McV in the MID model leads to 
an estimate of -477 McV < A'sat.2 < -241 MeV. 



While the estimated value of -477 MeV < A:sat,2 < 
—241 MeV in the present work has a small overlap with 
the constraint of Kt — —550 ± 100 MeV obtained in Ref. 

, [l9t from recent measurements of the isotopic depen- 
dence of the GMR in even-A Sn isotopes, the magnitude 
of the constrained Kr is still significantly larger than that 
of -477 MeV < A:sat,2 < -241 MeV. Recently, there are 
several works [4^,[53| on extracting the value of the A'sat,2 
parameter based on the idea initiated by Blaizot and col- 
laborators that the values of both ATo and Arsat.2 should 
be extracted from the same consistent theoretical model 
that successfully reproduces the experimental GMR en- 
ergies of a variety of nuclei. These studies show that 
there is no a single model (interaction) which can simul- 
taneously describe correctly the recent measurements of 
the isotopic dependence of the GMR in even-A Sn iso- 
topes and the GMR data of nuclei ^°Zr and "^^^Ph, which 
makes it difficult to accurately determine the value of 
Arsat.2 from the experimental GMR energies of a vari- 
ety of finite nuclei. As pointed out in [s^l, these features 
seem to sug gest that the K^ = — 550 ± 100 MeV obtained 
in Ref. [l8l . [19j may suffer from the same ambiguities al- 
ready encountered in earlier attempts [6l| to extract the 
ATo and Arsat.2 of infinite matter from finite-nuclei extrap- 
olations. This problem remains as an open challenge, and 
both experimental and theoretical insights are needed in 
the future. 
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